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ABSTRACT 



Context. The steep-spectrum radio quasars (SSRQs) are powerful radio sources, with thermal emission from accretion disk and jet 
nonthermal emission likely both contributing in the Ultraviolet (UV)/optical luminosity, however the former may play a dominant 
role. While the UV/optical variability of SSRQs has been poorly studied, little is known on the mechanism of their variability. 
Aims. We investigate the mechanism of the UV/optical variability of SSRQs. 

Methods. A sample of eighteen SSRQs has been established in SDSS Stripe 82 region in our previous works, in which the flux and 
spectral variability have been studied. In this work, we construct the flux-flux diagram using SDSS u and i multi-epoch data for these 
eighteen SSRQs. The standard accretion disk model is used to fit the flux-flux variations, in order to explore the variability mechanism. 
Results. The model fit to flux-flux diagram are tuned with fixed black hole mass and varying accretion rate. We found that the flux- 
flux diagram of all our SSRQs can be qualitatively described by the standard accretion disk model with change of the accretion rate. 
Although non-thermal jet power-law emission can also qualitatively reproduce the variability, the reasonable accretion rates and black 
hole masses required to fit the flux-flux variations suggest that a disk emission with variable accretion rate is a plausible description 
of the data. 
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1. Introduction 

Although radio loud Active galactic nuclei (AGNs) only rep- 
resent 10% of whole AGNs population (e.g. Kellermann et al. 
1989), they are important in AGNs study, for example, to study 
the jet production, composition, collimation and even the feed- 
back to the environment. As two subsets of radio-loud AGNs, 
flat-spectrum radio quasars (FSRQs) and steep-spectrum radio 
quasars (SSRQs) are different in many respects. FSRQs together 
with BL Lac objects, so called Blazars, are the most extreme 
class of active galactic nuclei (AGNs), characterized by strong 
and rapid variability, high polarization, and apparent superlumi- 
nal motion. These extreme properties are generally interpreted 
as a consequence of non-thermal emission from a relativistic 
jet oriented close to the line of sight, which can well explain 
the large core-dominance in FSRQs in general. In contrast, the 
SSRQs are usually lobe-dominated radio quasars, and the radio 
lobe emission dominates over the radio core emission. Their jets 
are viewed at larger angles than blazars, hence, the beaming ef- 
fects of jets are not severe (see e.g. Liu et al. 2006), and the jet 
emission is not expected to dominate at optical bands (e.g. Gu & 
Ai2011a,b). 

AGNs exhibit variability at almost all wavelengths (e.g. 
Wiita 1996; Vanden Berk et al. 2004). Multiwavelength studies 
of variations in the radiation emitted from AGNs have played im- 
portant roles in exploring the physical conditions near the cen- 
ter of AGNs. There have been extensive studies of the optical 
variability of blazars (e.g. Ghisellini et al. 1997; Fan et al. 1998; 
Massaro et al. 1998; Ghosh et al. 2000; Clements & Carini 2001 ; 
Raiteri et al. 2001; ViUata et al. 2002; Vagnetti et al. 2003; Wu 
et al 2005, 2007; Gu et al. 2006; Hu et al. 2006; Poon et al. 



2009; Rani et al. 2010; Gu & Ai 201 la). While it is generally ac- 
cepted that the nonthermal emission from relativistic jet oriented 
close to the line of sight dominates the optical continuum, the 
situation seems more complicated in FSRQs, with the evidence 
of thermal emission in FSRQs (e.g. Plan et al. 1999; Grandi & 
Palumbo 2004; Raiteri et al. 2007; D' Ammando et al. 201 1), and 
the redder-when-brighter trend found in several FSRQs (e.g. Gu 
et al. 2006; Rani et al. 2010; Gu & Ai 201 la; Wu et al. 201 1). 

Intrinsic ally, the variations in the brightnesses of AGNs are 
generally caused by physical variations in the jet and accretion 
disk. However, the contribution of each component to the ob- 
served variability varies from source to source. The properties 
of steep-spectrum radio quasars (SSRQs) are in-between those 
of both FSRQs and radio-quiet quasars. Therefore, it is not clear 
whether the jet nonthermal or accretion thermal emission is pri- 
marily responsible for the variation in SSRQs, even the situation 
can be more complicated if two components are mixed together 
with comparable contribution. However, the optical and spectral 
variations of SSRQs have been poorly studied, and few expla- 
nations have been presented (e.g. Stalin et al. 2004, 2005; Gu & 
Ai 2011b), in comparison to many investigations in radio-quiet 
AGNs (e.g. Stahn et al. 2004, 2005; Gupta & Joshi 2005; Wilhite 
et al. 2005; Ai et al. 2010) and FSRQs. The optical variability of 
SSRQs has been recently investigated by Gu & Ai (201 lb) for a 
sample of eighteen SSRQs assembled from SDSS Stripe 82 re- 
gion. As shown in Gu & Ai (201 lb), the Mg II line to continuum 
luminosity ratio and the anti-coiTelation between the variability 
at r band and the Eddington ratio of SSRQs are all similar to 
those of radio-quiet AGNs. These results strongly argued that 
the thermal emission from accretion disk might be the domi- 
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nant one in the optical continuum, and it may be responsible for 
the variability of SSRQs. In terms of the long-term variability, 
the change in accretion rate is used to explain the optical vari- 
ations in radio quiet AGNs (e.g. Li & Cao 2008). In this work, 
we investigate whether the change of accretion rate in standard 
accretion disk model can well describe the Ultraviolet/optical 
variabiHty of SSRQs. 

The layout of this paper is as follows: in Section 2, we de- 
scribe the source sample; the results of model fit to the flux varia- 
tions are outlined in Section 3; Section 4 includes the discussion; 
and in the last section, we draw our conclusions. The cosmolog- 
ical parameters Hq - 70 km s"' Mpc"' , Qm = 0.3, and Qa = 0.7 
are used throughout the paper, and the spectral index a is defined 
as /v oc v^", where fy is the flux density at frequency v. 



2. Sample 

Our sample of eighteen SSRQs consists of five SSRQs in Gu 
& Ai (2011a) and thirteen SSRQs in Gu & Ai (2011b). The 
initial quasar sample was selected as those quasars in both the 
SDSS DR7 quasar catalogue (Schneider et al. 2010) and Stripe 
82 region. The SDSS DR7 quasar catalogue consists of 105,783 
spectroscopically confirmed quasars with luminosities brighter 
than Mi - -22.0, with at least one emission line having a full 
width at half-maximum (FWHM) larger than 1000 km s"' and 
highly reliable redshifts. The sky coverage of the sample is about 
9380 deg^ and the redshifts range from 0.065 to 5.46. The five- 
band (m, g, r, /, z) magnitudes have typical errors of about 0.03 
mag. The spectra cover the wavelength range from 3800 A to 
9200 A with a resolution of ^ 2000 (see Schneider et al. 2010, 
for details). The Stripe 82 region extending over right ascension 
a = 20'' - A^ and declination -1°.25 < 6 < -h1°.25, was re- 
peatedly scanned during the SDSS-I phase (2000 - 2005) under 
generally photometric conditions. This region was also scanned 
repeatedly over the course of three 3-month campaigns in three 
successive years in 2005 - 2007, as part of the SDSS Supernova 
Survey (SN survey). The data are well-calibrated (Lupton et al. 
2002). 

We cross-correlate the initial quasar sample with the Faint 
Images of the Radio Sky at Twenty centimeters (FIRST) 1 .4- 
GHz radio catalogue (Becker, White & Helfand 1995), the Green 
Bank 6-cm (GB6) survey at 4.85 GHz radio catalogue (Gregory 
et al. 1996), and the Parkes-MIT-NRAO (PMN) radio continuum 
survey at 4.85 GHz (Griflith & Wright 1993). The initial quasar 
sample was first cross-correlated with the SDSS quasar positions 
and the FIRST catalogue to within 2 arcsec (see e.g. Ivezic et al. 
2002; Lu et al. 2007). The resulting sample of SDSS quasar po- 
sitions was then cross-correlated with both the GB6 and PMN 
equatorial catalogues to within 1 arcmin (e.g. Kimball & Ivezic 
2008). Owing to the different spatial resolutions of FIRST, GB6, 
and PMN, single or multiple FIRST counterparts were found to 
within 1 arcmin, although only a single GB6 and/or PMN coun- 
terpart existed. The radio spectral index a^ was then calculated 
between the single or integrated FIRST and/or NRAO VLA Sky 
Survey (NVSS) 1.4 GHz within 1 arcmin and either or both of 
the GB6 and PMN 4.85 GHz (see details in Gu & Ai 201 la,b). 
We define a quasar to be a SSRQ according to its radio spectral 
index a, > 0.5. Finally, five SSRQs with single FIRST counter- 
part to within 1 arcmin of the SDSS positions were presented in 
Gu & Ai (2011a), while thirteen SSRQs with multiple FIRST 
counterparts were listed in Gu & Ai (201 lb). To have the radio 
spectral information, we decided to consider only the sources 
with a radio counterpart in both 1.4 and 4.85 GHz radio images. 



In spite of this, our final sample, although not complete, still in- 
cludes a quite large number (18) of SSRQs. 

The source sample is listed in Table [1] in which the black 
hole mass and the Eddington ratio from Gu & Ai (2011b) are 
shown. In calculating the Eddington ratio, the bolometric lumi- 
nosity is estimated as LboI = IOLblr (Netzer 1990), in which 
the broad line region (BLR) luminosity Lblr is derived follow- 
ing Celotti et al. (1997) by scaling the strong broad emission 
lines Hy6, Mg II, and C IV to the quasar template spectrum of 
Francis et al. (1991), in which Lya is used as a flux reference of 
100. 

The multi-epoch photometric observations in the Stripe 82 
region enable us to investigate the optical variability of the se- 
lected quasars. As in Gu & Ai (2011a, b), we directly used the 
point-spread-function magnitudes in the Catalog Archive Server 
(CAS) Stripe82 database from the photometric data obtained 
during the SDSS-I phase from data release 7 (DR7; Abazajian 
et al. 2009) and the SN survey during 2005 - 2007. We selected 
the sources classified as point sources in all observational runs, 
and only data with good measurements (high-quality photome- 
try) were selected, of which the ugriz magnitude are required to 
be brighter than the magnitude limit. Finally, the data taken in 
cloudy conditions were also excluded. 



3. Results 

Besides the conventional flux versus color relation, the flux-flux 
relation has been applied to examine the color variability for 
quasars (e.g. Sakata et al. 2011; Schmidt et al. 2012). In this 
work, to obtain a wide wavelength coverage in investigating the 
spectral variability, the u and / magnitude in the same night are 
used to plot flux-flux diagram. The reason of using / magnitude 
instead of the z magnitude as the longest wavelength is that the 
signal-to-noise ratio is much higher in the former than that of the 
latter for most targets. The corresponding wavelength of SDSS 
M and / waveband in the source rest frame are shown in Table [T] 
The Eddington ratio LBoi/^Edd of all SSRQs is -2.0 < 
log LBoi/ifidd < except for three sources (SDSS J213004.75- 
010244.4, SDSS J012517.14-001828.9 and SDSS J015832.51- 
004238.2, see Table [1]), which is indicative of a standard thin 
disk (Shakura & Sunyaev 1973) being present in SSRQs. The 
low LBoi/^Edd in SDSS J213004.75-010244.4, and high values in 
SDSS J012517.14-001828.9 and SDSS J015832.51-004238.2, 
make the standard thin disc hardly applicable. However, the esti- 
mation of Eddington ratio is subject to large uncertainties both in 
black hole mass and BLR luminosity (see section 4). The stan- 
dard accretion disk model used in Li & Cao (2008) has success- 
fully explained the correlation between the optical-UV variabil- 
ity amplitude and black hole mass found in Wold et al. (2007) 
and Wilhite et al. (2008). In this work, we directly adopted the 
model of Li & Cao (2008) except that the temperature distribu- 
tion (the Equation (1) in Li & Cao 2008) was modified as 






(1) 



(Frank, King & Raine 2002). The temperature distribution in Li 
& Cao (2008) represents the disk temperature on the equatorial 
plane. In contrast, the temperature distribution we used here (in 
Equation (1)) describes the radial dependence of the effective 
temperature derived from the emergent flux from an accretion 
disk, which is more appropriate for our purpose. We calculated 
accretion-disk spectra assuming a steady geometrically thin, op- 
tically thick accretion disk (Shakura & Sunyaev 1973). In this 
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case the emitted flux is independent of viscosity, and each el- 
ement of the disk face radiates roughly as a blackbody with 
a characteristic temperature that depends only on the mass of 
the black hole, Mbh, the mass accretion rate, M, and the ra- 
dius of the innermost stable orbit (Peterson 1997; Frank, King 
& Raine 2002). We have adopted the Schwarzschild geometry 
(non-rotating black hole), and for this the innermost stable or- 
bit is at Tin = 6rg, where r^ is the gravitational radius defined 
as Tg - GMbh/c^, G is the gravitational constant, and c is the 
speed of light. Unlike blazars, SSRQs are expected to have the 
jet pointed away from the observer We have assumed that the 
disk is viewed with an inclination angle of 20°, which is the 
mean viewing angle of lobe-dominated radio quasars in Gu, Cao 
& Jiang (2009). The accretion-disk spectrum is fully constrained 
by the two quantities, the mass accretion rate and the black hole 
mass. 

We fit the accretion disk model to the flux-flux plots of in- 
dividual SSRQs with a free parameter of black hole mass and 
by changing the mass accretion rate M. The flux-flux diagrams 
are shown in Fig. [T]represented by the rest frame flux density at 
frequencies corresponding to /l„ and /I,, with dotted lines repre- 
senting the best-fit accretion disk model. The best-fit black hole 
mass, reduced x^ values and the range of the accretion rate m 
required to reproduce the flux-flux variations are listed in Table 
[T] The accretion rate m is defined as m - M/Mem, in which M is 
the mass accretion rate from model fit, and MEdd is the Eddington 
accretion rate (the Eddington luminosity Lndd - O.lMEddC^). 
Assuming Lboi = O.lMc^, the accretion rate m is equal to 
the Eddington ratio Lboi/i^Edd- The standard thin disk seems to 
roughly describe the observed flux to flux variations. However, 
the;t'^-test values are generally large. 

As SSRQs are sources with a jet component, we also tested 
the possibility that a non-thermal power law emission can cause 
the observed flux-flux variations. As in the jet dominated blazars, 
the bluer-when-brighter trend is commonly observed, a power- 
law function of /„ = axf was adopted (see Sakata et al. 201 1). 
The fitting of the power-law function to the flux-flux variations 
are presented in Fig. [T] and the best-fit parameters and the re- 
duced ;t'^ values are listed in Table [T] We found that the reduced 
X^ value of the power-law fit is comparable to that of the standard 
thin disk for all quasars, although slightly improved. It seems 
difficult to distinguish the thermal disk and nonthermal jet emis- 
sion from a statistical point of view. 



4. Discussions 

While FSRQs are usually associated with core-dominated ra- 
dio quasars, SSRQs are generally related to lobe-dominated 
ones, usually with two large-scale optically thin radio lobes. The 
beaming effect is usually not strong in SSRQs because of the rel- 
atively large viewing angle. As shown in Gu & Ai (2011b), the 
Mg II line to continuum luminosity ratio and the anti-correlation 
between the variability at r band and the Eddington ratio of 
SSRQs are all similar to those of radio quiet AGNs. These results 
strongly argued that the thermal emission from accretion disk 
might be the dominant one in the optical continuum, and it may 
be responsible for the variability of SSRQs. While the change 
of accretion rate can qualitatively explain the optical variations 
in radio quiet AGNs (e.g. Li & Cao 2008), it is expected that 
it can also be applied to SSRQs. From a qualitatively point of 
view, an accretion rate variations seems to work also for SSRQs. 
However, according to the ;t'^-test, only two sources are consis- 
tent with this expectation. 



Quasars usually have strong broad emission lines (BELs), 
which could affect the photometry when they are in the SDSS 
u and / bands. Following Elvis et al. (2012), we found that only 
three strongest BELs have contribution larger than 3%, i. e. Lya, 
Hff and li/3, for a typical bandwidth of 1000 A for u and / bands 
(see equation (1) and Table 2 of Elvis et al. (2012)). Therefore, 
we only consider these three lines. In our sample, four sources 
(labelled in Table 1) are contaminated by BELs, of which Lya 
affects in two sources (SDSS J012401.76-H003500.9 and SDSS 
J213513. 10-005243. 8), Ha in SDSS J235156.12-010913.3, 
and Hj6 in SDSS J021225.56-I-010056.1, respectively. In princi- 
ple, the line contribution to the u and / fluxes could be subtracted 
considering the SDSS spectra or the method proposed by Elvis et 
al. (2012). However, we prefer to not remove any Une contribu- 
tion. It is quite improbable that the large discrepancy observed 
between data and standard thin disc model can be entirely as- 
cribed to emission line features. Note that a large reduced x^ 
(6.4) is also reported for SDSS J022508.07-I- 001707.2, a SSRQ 
in the sample without BELs contaminations. 

Since the reduced ^^ values are generally large, it is possible 
that the accretion disk model is oversimplified and that variations 
of the accretion rate are not sufficient to explain the observed flux 
variability. As argued by Zuo et al. (2012), the change of accre- 
tion rate is important for producing the observed optical variabil- 
ity, but other physical mechanisms still need to be considered in 
modifying the simple accretion model for quasars. Schmidt et 
al. (2012) analyzed a large sample of 9093 quasars from SDSS 
Stripe 82 region, and claimed that on timescales of years quasar 
variability does not reflect changes in the mean accretion rate 
and some other mechanism must be at work, presumably some 
disk instability. The authors argued that this picture was con- 
firmed by the comparison of the observed color variability to 
sequences of steady-state accretion disk models with varying ac- 
cretion rates, which exhibit much less color variability as a func- 
tion of accretion rate. While this is claimed for quasar ensem- 
bles, the thin accretion disk models could well match the data 
both in color and in the change of color with changing luminos- 
ity by changing the accretion rate in some individual objects (see 
Fig. 11 in Schmidt et al. 2012). In contrast, Sakata et al. (2011) 
found that the multi-epoch flux-to-flux plots could be fitted well 
with the standard accretion disk model changing the mass ac- 
cretion rate with a constant black hole mass for luminous radio- 
quiet quasars. However, their sample size is rather small, only 
containing ten sources. Moreover, Pereyra et al. (2006) found 
that the composite differential spectrum of two epochs of ob- 
servations for hundreds of SDSS QSOs in the rest-frame wave- 
length range 1300 - 6000 A can be successfully fitted by the 
standard thin disk model, provided that their accretion rates vary 
from one epoch to the next. It shoud be noted that the composite 
spectrum only represents the average characteristics of QSOs. It 
thus is unclear whether the scenario is applicable to the individ- 
ual objects. Our results, although based on a smaller Stripe 82 
region sample are in agreement with Schmidt et al. (2012). 

As our sources are radio emitters, we can not exclude that the 
jet contributes to the optical/UV emission. The large reduced ;if^ 
values then may naturally reflect the contribution of jet emission 
in the UV/optical band, in addition to the thermal accretion disk 
emission. The nonthermal contribution to the MIR luminosity 
has been investigated for a sample of powerful 3C RR galax- 
ies and quasars with the IRS and MIPS instruments on Spitzer 
(Cleary et al. 2007). Due to the low radio frequency selection, 
their quasars are mostly steep spectrum sources, including sev- 
eral compact steep spectrum sources (CSS). By fitting the con- 
tinuum with synchrotron (jet and lobe) and dust components. 
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they found that nonthermal processes can contribute a signifi- 
cant proportion (up to 90%) of the measured infrared emission in 
some quasars, although they also shown that the nonthermal con- 
tribution may not be severe in the sources with lower core dom- 
inance (see their Fig. 7), which is typical in SSRQs. However, 
our disk model results (see Table 1) indicate that the accretion 
rates of the studied sources are in general agreement with values 
expected for standard accretion disks. 

According to the Eddington ratio estimated using the bolo- 
metric luminosity from BLR, two sources have Eddington 
ratio larger than 10*^^ ~ 3, and one source SDSS 
J2 13004.75 -01 0244.4 has the Eddington ratio smaller than 0.01. 
On the contrary, the disk model fitted accretion rate of these 
three sources are all well in the range of standard thin disk model 
> 0.01 . The m difference is likely due to the large uncertainties in 
the black hole mass estimated from BLR luminosity and the em- 
pirical relation. We estiamte the errors of the model fit black hole 
masses from the statistical uncertainty of ^^ fitting of the model, 
which is listed in Table [1] While the error on the Mbh estimates 
from BLR is ~ 0.4 dex (e.g. Vestergaard & Peterson 2006), the 
uncertainties of the black hole mass from accretion disk model 
fit is much less. The black hole masses estimated from two meth- 
ods are compared in Fig. |2]f or 17 SSRQs, after excluding the low 
accretion rate source SDSS J013514.39-000703.8. Generally, a 
good consistence is found between these two masses. Therefore, 
in spite of the fact that the disk model fits are not acceptable for 
most of the sources and a power law can also reproduce the ob- 
servations, the reasonable accretion rates and black hole masses 
required to fit the flux-flux variations strongly suggest that a disk 
emission with variable accretion rate is a plausible description of 
the data. 

Optical polarimetry could be extremely useful in investigat- 
ing the nature of optical/UV emission in radio sources. As shown 
by the recent studies of 3C 120 and 3C 111, the low optical 
polarization supports the thermal origin of the optical emission 
(Chatteijee et al. 2009, 2011), however, the detection of polar- 
ization percentages of 3% or even higher in 3C 111 suggests 
that during particular activity phases an additional contribution 
from non-thermal emission from the jet may be present (Jorstad 
et al. 2007; Chatterjee et al. 2011). Motivated by this fact, the 
polarimetry observation are therefore needed to further explore 
the jet contribution in the optical band for SSRQs, in addition to 
the photometric and spectroscopic monitoring. 

5. Summary 

We explore the flux-flux diagram using SDSS u and / multi- 
epoch data for a sample of eighteen SSRQs in SDSS Stripe 82 
region established in our previous works. The standard accre- 
tion disk model is used to fit the flux-flux variations, in order 
to explore the variability mechanism. We found that the flux- 
flux diagram of all our SSRQs can be qualitatively described 
by the standard accretion disk model with change of the accre- 
tion rate. Although non-thermal jet power-law emission can also 
qualitatively reproduce the variability, the reasonable accretion 
rates and black hole masses required to fit the flux-flux varia- 
tions suggest that a disk emission with variable accretion rate is 
a plausible description of the data. 
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Table 1. Source list: Col. 1 - SDSS source name; Col. 2 - redshift; Cols. 3 - 4 - the corresponding wavelength of SDSS u and / wavebands in the source rest frame, respectively 
(e, /, g jjjg prominent broad emission lines Lya, H/3, and Ha in the photometric wavebands, respectively); Col. 5 - " black hole mass in unit of solar mass estimated from the 
empirical relation of BLR radius and the BLR line width; Col. 6 - the Eddington ratio I - L-^oi^jL^uD with the black hole mass in Col. 5, and the bolometric luminosity estimated 
from BLR luminosity; Cols. 7 - 8 - * black hole mass in unit of solar mass, and the range of the accretion rate from the model fit with a standard accretion disk model, respectively; 
Col. 9 - '^ the reduced ^^ of the accretion disk model fit; Col. 10 -'' the reduced ;^^ of the power-law fit; Cols. 11 - 12 - the parameters of the power-law fit /„ -axf . 
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Fig. 1. The flux-flux diagram for our SSRQs represented by the 
rest frame flux density at frequencies corresponding to /l„ and /I, 
(see Table[T]l fory- and jc-axis, respectively. The dotted Hne is the 
best-fit standard accretion disk model with varied mass accretion 
rates and a constant black hole mass, and the dash-dotted line 
represents the best-fit of a power-law. The solid circles with error 
bars represent the multi-epoch flux data. 
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Fig. 1. - continued. 
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Fig. 1. - continued. 
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Fig. 2. The black hole mass estimated from the empirical 
relation of BLR radius and the BLR line width compared 
with that of the best-fit standard accretion disk model for 17 
SSRQs after excluding the source with low accretion rate SDSS 
J013514. 39-000703. 8. The solid triangles are four sources with 
strong emission line contamination in either u or ; bands. 



10.5 



